2009 Fermi Symposium, Washington, B.C., Nov. 2-5 



Initial Results of New Tomographic Imaging of the Gamma-Ray Sky 
with BATSE 

G. L. Case, M. L. Cherry 

Dept. of Physics and Astronomy, Louisiana State University Baton Rouge, LA 70803, USA 
J. C. Ling, M. Lo, T. Shimizu 

Jet Propulsion Laboratory, California Inst, of Technology, Pasadena, CA 91109, USA 
Wm. A. Wheaton 

Infrared Processing and Analysis Center, California Inst, of Technology, Pasadena, CA 91 125, USA 



We describe an improved method of mapping the gamma-ray sky by applying the Linear Radon Transform 
to data from BATSE on NASA's CGRO. Based on a method similar to that used in medical imaging, we use 
the relatively sharp (~ 0.25°) limb of the Earth to collimate BATSE's eight Large Area Detectors (LADs). 
Coupling this to the ~ 51-day precession cycle of the CGRO orbit, we can complete a full survey of the sky, 
localizing point sources to < 1° accuracy. This technique also uses a physical model for removing many sources 
of gamma-ray background, which allows us to image strong gamma-ray sources such as the Crab up to ~ 2 
MeV with only a single precession cycle. We present the concept of the Radon Transform technique as applied 
to the BATSE data for imaging the gamma-ray sky and show sample images in three broad energy bands 
(23-98 keV, 98-230 keV, and 230-595 keV) centered on the positions of selected sources from the catalog of 130 
known sources used in our Enhanced BATSE Occultation Package (EBOP) analysis system. Any new sources 
discovered during the sky survey will be added to the input catalog for EBOP allowing daily light curves and 
spectra to be generated. We also discuss the adaptation of tomographic imaging to the Fermi GBM occultation 
project. 



1. Introduction 

The soft gamma-ray sky in the 200 — 2000 keV range 
is one of the more poorly studied regions of the elec- 
tromagnetic spectrum. Many of the most intrinsically 
interesting objects in the universe emit much of their 
primary energy in this regime. The last complete all- 
sky survey above ~ 100 keV was made by HE AO in 
the late 1970s and early 1980s. While missions such as 
INTEGRAL have considerably better angular resolu- 
tion and sensitivity than HE AO A4, they have a small 
field of view and reduced sensitivity above 200 keV. 
Because of the penetrating power of gamma rays, mas- 
sive collimators would be needed in order to achieve 
reasonable sensitivity in the ~ 200 — 2000 keV region. 




Figure 1: Left: Time history from a BATSE LAD. 'A' 
marks the time at which the Crab is occulted, 'B' marks 
the time during which the Crab is blocked by the Earth, 
and 'C marks the time the Crab reappears. Right: 
Schematic of the Earth occultation process. 



The Burst And Transient Source Experiment 
(BATSE) instrument was one of the four instruments 
on the Compton Gamma Ray Observatory (CGRO). 
It consisted of 8 Nal large area detectors (LADs), 50 
cm in diameter and 1.27 cm thick, each paired with 
a Spectroscopy Detector, 12.7 cm in diameter and 
7.62 cm thick, located on the corners of the space- 
craft. Although the LADs have no direct imaging ca- 
pability, it is possible to monitor known gamma-ray 
sources using the Earth occultation technique [TJ [2] . 
The Earth occultation technique utilizes the relatively 
sharp (~ 0.25°) limb of the Earth to collimate photons 
from individual sources. Figure [I] demonstrates how 
the occultation technique works. The count rate mea- 
sured by a LAD is comprised of contributions from 
all of the cosmic sources not occulted by the Earth 
as well as a background component. When a source 
is occulted by the Earth, the count rate measured in 
the LAD drops, producing a downward step-like fea- 
ture in the count rate. When the source reappears, an 
upward step is produced. 

The Enhanced BATSE Occultation Package 
(EBOP) was developed at JPL pQ to fit the con- 
tinuous count rate (CONT) data for an entire day 
to a model including all of the occultation steps 
that occur for the sources in the input catalog along 
with a semi-physical background model. The output 
consists of the average daily count rates in each LAD 
for each source in each of 16 energy bands. However, 
EBOP requires a priori knowledge of the positions 
of the gamma-ray sources, i.e. a predetermined 
input catalog. The discovery of new sources must be 
accomplished using other methods, such as imaging. 

Using BATSE, images of the sky in the 25 - 2000 
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keV energy range can be obtained utilizing tomo- 
graphic imaging techniques such as the Linear Radon 
Transform, a technique used in medical imaging ap- 
plications. Applying this technique to the BATSE 
CONT data using the Earth limbs gives an angular 
resolution roughly equal to the ~ 15' size of the atmo- 
sphere, as viewed from CGRO, that the source must 
traverse as it occults. With nine years of data avail- 
able, images of the entire sky can be generated with 
high sensitivity, and the BATSE dataset can be used 
to search for previously unknown gamma-ray sources. 



2. The Linear Radon Transform 

The Linear Radon Transform (LRT) is the integral 
transform of an image / across all lines L in the (x, y) 
plane into the (p, 0) plane, where 9 is the angle that 
the line L makes with respect to the x-axis and p is 
the perpendicular distance from the line to the origin, 

R(I) = J I(x,y) dL = F(6,p), (1) 

where L = {y = x tan + pj cos 0}. 

In Fig. [2) the LRT and its inverse are demonstrated 
for a sample image. The middle image shows how the 
features in the (x, y) plane are transformed into curves 
in the (p, 0) plane. The Inverse LRT transforms the 
integrated curves back into the (x, y) plane. 

This technique was first applied to Earth occupa- 
tion using BATSE data below 100 keV by Zhang et 
al. [3J. We implement it as follows. The 2.048-second 
CONT data is rebinned to 16.384 seconds and the 
background model, as obtained with EBOP, is sub- 
tracted. The difference in the counts measured by 
a detector in two adjacent time bins, to and t\ with 
ti > to for rising steps and t\ < to for setting steps, 
gives the integrated number of counts from the arc 
defined by the position of the Earth's limb at time 
to. Thus an Earth occultation measurement is a LRT 
of the sky across a single "line" . The LRT assumes a 
straight line, so the image tile must be restricted to 
small enough regions of the sky that the Earth's limb 
can be approximated as a straight line. As the plane of 



the CGRO orbit precesses (with a period of ~ 51 days 
given its 28.5° inclination), the Earth's limb will be 
oriented at different angles with respect to the image 
tile. An image is produced by performing the Inverse 
LRT on all of the limbs in one precession period. 

3. The Earth Limbs and the Point Spread 
Function 

The point spread function (PSF), and hence the ex- 
pected shape of the image, depends on the declination 
of the source. Figure [3] shows the limb orientations 
(drawn as multi-colored lines) that are sampled dur- 
ing a complete orbit plane precession cycle for sources 
at three different declinations, along with the simu- 
lated PSFs and the actual images. For a source at 
low declination (e.g. 3C273), fewer limb orientations 
are sampled (that is, the lines L have a smaller range 
of angles, 0) resulting in a PSF that is elongated in 
the north-south direction. For sources at modest dec- 
lination (e.g. IE 1740-29), about half of the possible 
limb orientations are sampled, resulting in a PSF that 
is slightly elongated in the north-south direction with 
an 'x' pattern in the wings. For sources with declina- 
tion \S\ > 45° (e.g. GX 301-2), nearly all of the limb 
orientations are sampled and the PSF is roughly circu- 
lar or even elongated in the east-west direction. The 
PSFs are calculated by simulating the response to a 
unit intensity source for a particular time period, tak- 
ing into account the expected limb orientations and 
the gaps that occur in the real data. Gaps can ap- 
pear in the data due to events such as passes through 
the South Atlantic Anomaly (when the detector high 
voltage is turned off), gamma-ray bursts, transmission 
dropouts, etc. As can be seen in Fig. |3j the shapes of 
the images agree well with the simulated PSFs. 



4. Images 

Figure [4] shows a selection of images generated with 
the Inverse LRT in the 23-98 keV (band 1) energy 
band, while Fig. [5] and Fig. [6] show images from the 
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Figure 3: Top: Orientation of limbs during one precession cycle for three sources at different declinations. The colors of 
the lines have no meaning and are used to distinguish individual lines. Middle: Simulated PSFs for the same three 
sources. Bottom: Actual images in the 23-98 keV band of the three sources averaged over 20 precession cycles. 



same time period for the 98-230 keV (band 2) and 
230-595 keV (band 3) energy bands, respectively. The 
images for persistent sources are summed over 20 pre- 
cession periods, from TJD 8393-9397, while transient 
sources are summed over one or more complete preces- 
sion cycles spanning the time of the outburst. Com- 
plete precession cycles are used in order to make sure 
that all available limb orientations are sampled, which 
creates a more uniform image. All detectors that view 
the source at an angle <p < 60° from the detector nor- 
mal are combined, weighted by the detector efficiency. 

The Crab pulsar/nebula is the brightest source seen 
by BATSE in all energy bands. It is a persistent source 
and is clearly seen in all three energy bands, even over 
a single precession cycle. In fact, the Crab is also seen 
in band 4 (595-1800 keV) gj. Cyg X-l, a high-mass x- 
ray binary harboring a black hole, is a persistent but 
variable source that is also seen clearly in all three en- 
ergy bands. GRO J0422+32 is a transient source pow- 
ered by a black hole that was discovered by BATSE, 
and at its peak was several times brighter than the 



Crab in band 1. During its outburst, a hard spec- 
tral component was observed by Ling & Wheaton [5] , 
which is evident in the band 3 image. 

GRO J1719-24 is another transient black hole 
source discovered by BATSE. During its outburst, it 
also was observed to have a hard spectral component 
[6] and is about the same intensity as the Crab in band 
2. GX 339-4 is also a transient black hole source that 
underwent four outbursts from 1991 to 1994. The GX 
339-4 images in Fig. [4] and Fig. [5] are the averages 
of the four individual images during the individual 
outbursts. The band 2 image shows clear evidence 
of emission at the position of GX 339-4. The black 
hole transient GS 1354-654 underwent an outburst in 
1997-98. Though it was not as bright as the other 
black hole transients, the image of GS 1354-645 shows 
a significant peak at the known position of the source 
(as marked by the sign). 

Sco X-l and Vela X-l are low-mass x-ray bina- 
ries (LMXBs) that are persistent but highly variable 
sources with soft gamma-ray spectra. They are mod- 
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Figure 4: Images in the 23-98 keV band (band 1). The intensity units are arbitrary and the ' + ' sign marks the known 
positions of the sources. 



erately bright in band 1, but neither of these sources 
are seen definitively at energies above 100 keV [4J. 
GRS 1915+105 and 4U 1700-37 are both LMXBs 
(GRS 1915+105 containing a black hole and 4U 1700- 
37 containing a neutron star). Both sources are mod- 
erately bright in band 1 but much fainter in band 2, 
and neither source is seen above 230 keV [4]. GRS 
1009-45 is a transient black hole source that under- 
went a single outburst in 1993. At its peak, it was 
nearly as bright as the Crab in band 1, but much 
fainter in band 2, and not visible above 230 keV [4]. 
The extragalactic source NGC 4151 is a relatively faint 



source that shows up consistently in band 1 images 
summed over fewer precession cycles, with the cen- 
tral feature of the image consistent with the known 
position. 

The shapes of the features in the images are con- 
sistent with expected PSFs based on the source decli- 
nation as discussed in Section 03 The low declination 
sources (GRS 1915+105 and Sco X-1) are very elon- 
gated in the north-south direction, while the moder- 
ate declination sources (Crab, GRO J0422+32, GRO 
J1719-24, Cyg X-1, and 4U 1700-37) are only sightly 
elongated with the characteristic 'x' shape further 



eConf C091122 



2009 Fermi Symposium, Washington, B.C., Nov. 2-5 



5 



24. 



~23| 

~22| 

o 

LU 

°21| 



Crab Nebula 



85 84 83 82 
R.A. (Deg) 



Cygnus X-1 



GROJ0422+32 



J0.5 




301 300 299 298 
R.A. (Deg) 




67 66 65 64 
R.A. (Deg) 



-23. 



~24| 

m 

Q 

~-25| 

U 

LU 

° -261 



GROJ1719-24 



I 



GX339-4 



262 



260 
R.A. (Deg) 



258 



1.5 

0.5 





258 256 254 
R.A. (Deg) 



-63 1 



a>-64| 
D 

d -65 1 

UJ 

a 
-661 



GS1354-645 



214 212 210 208 206 
R.A. (Deg) 



0.6 
0.4 
0.2 




Figure 5: Images in the 98-230 keV band (band 2). The intensity units are arbitrary and the '+' sign marks the known 
positions of the sources. 
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Figure 6: Images in the 230-595 keV band (band 3). The intensity units are arbitrary and the '+' sign marks the 
known positions of the sources. 



away from the peak. The high declination sources 
(GX 339-4, GS 1354-645, Vela X-1, and GRS 1009- 
45) are nearly circular or even elongated in the east- 
west direction. Weaker sources such as NGC 4151 and 
GS 1354-645 have shapes that are more heavily influ- 
enced by source confusion (contributions from a limb 
that has another bright source occulting at the same 
time), other transient events in the data, e.g. gamma- 
ray bursts, solar flares, etc., or gaps in the data. 



5. Conclusions and Future Work 

We have briefly described the concept of Earth oc- 
cupation and the application of the Linear Radon 
Transform technique to the BATSE LAD data for 
imaging the gamma-ray sky. We have shown prelimi- 
nary images of cosmic gamma-ray sources at energies 



up to ~ 600 keV to demonstrate the feasibility and 
capability of the technique for generating images of 
the sky in the elusive soft gamma-ray regime, localiz- 
ing individual sources to < 1° accuracy. The images 
are consistent with the simulated PSFs showing the 
dependence of the shape of the features on the decli- 
nation. 

Ongoing work includes quantifying the significance 
of the features in the images and defining the sen- 
sitivity limits for this imaging method. We are also 
developing image deconvolution techniques in order to 
improve the quality of the images and more accurately 
locate the positions of individual sources, both those 
that are known and those that are newly discovered. 
Our long term goal is to create a map of the entire 
sky in the four broad energy bands 23-98 keV, 98- 
230 keV, 230-595 keV, and 595-1800 keV. This would 
be the first all-sky map in the ^ 200 — 2000 keV re- 
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gion since HE AO A4. The purpose of the sky map 
is to look for new sources that were not previously 
detected, which would then be added to our EBOP 
input catalog. The EBOP analysis system will then 
be rerun to both obtain fluxes for the new sources and 
to get more accurate fluxes for the previously known 
sources by reducing the systematic errors currently 
present due to the unknown sources in the sky. 

This work is directly applicable to the Gamma- 
Ray Burst Monitor (GBM) on Fermi, an instrument 
similar to BATSE. GBM can also be used to per- 
form Earth occultation measurements, and a program 
is already underway to monitor gamma-ray sources 
[7]. The imaging technique presented here using the 
LRT can be adapted to use with GBM, allowing new 
sources to be discovered. These new sources would 
be added to the input source catalog for continued 
monitoring, and the community would be alerted so 
that follow-up observations could be made at other 
wavelengths. 

The current method is based on the Linear Radon 
Transform which requires the image tiles to be rela- 
tively small (< 10°) in order to keep the "lines" of in- 
tegration (i.e. the Earth limbs) approximately linear. 
The geometry of the Earth occultation technique may 
be better described as an integration over the sphere 
with a spherical cap blocked by the Earth. The re- 
sulting integral is the spherical cap Radon transform 
which has many interesting properties, including the 
ability to create much larger image tiles. By utilizing 
this new transform in the imaging process, additional 
information may be extracted from the gamma-ray 
data to provide sharper images with more accurate 
flux estimates. Ultimately, we expect this enhanced 
technique to help identify new faint sources in the 
BATSE database and to be useful for current missions 
such as Fermi as well as future gamma-ray astronomy 
missions. 



The gamma-ray sky in the 200 keV to 10 MeV re- 
gion has yet to be fully explored to date. This is due 
in part to the complexity and cost of building an ap- 
propriate instrument, with wide field of view and good 
angular resolution for localizing sources to a fraction 
of a degree. The combination of a simple instrument, 
such as BATSE and GBM, with large omni-directional 
detectors and effective optics provided by the Earth 
limbs (or moon limbs, if flown in a lunar orbit), cou- 
pled to the application of tomographic imaging tech- 
niques, such as the LRT as shown in this work, would 
be an attractive, simple and cost effective option for 
future gamma-ray missions in the 200 keV to 10 MeV 
energy region. 
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